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^ ; Abstract 

Twisted mass lattice QCD (tmQCD), generalised to four Wilson quark flavours, can be 
■ used for the computation of some weak matrix elements related to A/ = 1/2 transitions. 

. Besides eliminating unphysical zero modes, tmQCD may alleviate the four-quark oper- 

^ \ ator renormalisation problems encountered in the calculation of CP-conserving if — > vr 

5h \ matrix elements with traditional Wilson fermion regularisation. With an active charm 

■ - - ' quark, the renormalisation of the K ^ it matrix elements requires at most the sub- 

traction of a linearly divergent counterterm. Furthermore, in the (partially) quenched 
approximation the twist angles can be chosen so that only a finite counterterm needs 
to be subtracted. 
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1 Introduction 



In the phenomenology of non-leptonic kaon and hyperon decays, the observed enhance- 
ment of decays with isospin change A/ = 1/2 over those with AI = 3/2 is referred to 
as the AI = 1/2 rule. For example, decays into a pion pair of definite isospin / can 
be parameterised through 



where 5j are the 5- wave vrvr phase shifts. The AI = 1/2 rule manifests itself in the 
large value of the ratio 



Although the exact origin of this enhancement is unclear, it is believed that it must be 
related to long-distance, non-perturbative physics. The contribution coming from en- 
ergy scales where perturbation theory can be reliably applied supplies only a small frac- 
tion of this ratio^. The rest should be an effect characteristic of low-energy QCD scales. 
Analytic studies have provided only a partial understanding of this non-perturbative 
enhancement factor pP . It is fair to say that we are still far from a quantitative predic- 
tion based on QCD. 

In principle, numerical simulations of lattice QCD are an ideal tool for the compu- 
tation of the required hadronic matrix elements. Several lattice regularisations of the 
fermionic action has been used to this aim. We cite indicatively recent work performed 
with Wilson [2], staggered [H], and domain wall fermions; for a more complete and 
fairly updated list of references, see ref. 0. In spite of the considerable effort invested 
in these simulations and several encouraging results, the ultimate goal is still far from 
being achieved, also because a number of conceptual and technical problems must be 
addressed beforehand. 

First, the finite spatial volume, characteristic of lattice simulations, strongly affects 
the two-pion final state [0]. In spite of this, it has been shown in [J] that physical K 
TTTT matrix elements can nevertheless be obtained from Euclidean correlation functions, 
computed in a finite volume (see also [H] for a related approach) . In particular a relation 
between the absolute values of the finite and infinite volume amplitudes is established. 
This, together with a separate measurement of the two-pion energies in finite volume, 
allows to fully reconstruct the desired decay amplitudes. However, the required spatial 
sizes L still tend to be larger by a factor of 2-3 than the volumes typically used for 
the numerical determination of the hadronic spectrum.^ Furthermore, unitarity is 
essential for the derivation of these results, so that the quenched approximation may 

^The precise value depends on the renormalisation scheme used in the computation; e.g. in typical 
MS schemes it roughly amounts to a factor of two. 

^For a recent proposal, mitigating this requirement, see ref. 




(1.1) 



Ao/^2| -22. 



(1.2) 
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not be applicable. Therefore, the practical feasibility of this programme hinges on 
possible theoretical developments, as well as further progress in simulation algorithms 
and computer hardware. 

The difficulty of computing the physical amplitudes directly has led to less am- 
bitious approaches, most of which aim at the determination of the coupling constants 
in the effective chiral Lagrangian TI?. This may be achieved by studying unphysical 
processes where no problem with final state interactions arises, such as K ^ n transi- 
tions, or K ^ vrvr decays with both pions at rest |in|lH12] . Once the effective coupling 
constants are known, the physical amplitudes can be inferred, up to higher order cor- 
rections in the chiral expansion. However, one might expect the next-to-leading order 
corrections to be large, especially if the matching to data from numerical simulations is 
performed at unphysically large meson masses. Assuming that the chiral expansion re- 
mains valid in this regime, it is thus desirable to include higher order terms in the chiral 
Lagrangian |13ll4j . Unfortunately, since this introduces many more unknown parame- 
ters already at next-to-leading order, the whole approach becomes quite cumbersome 
if not impractical. 

Whatever the chosen strategy, a difficult problem consists in the renormalisation 
of the four-quark operators which appear in the effective weak Hamiltonian. The de- 
tails strongly depend on the symmetries of the lattice regularisation; regularisations 
which preserve chiral symmetry are clearly advantageous |15piHpi7| . and are increas- 
ingly used despite their high computational costs. If computationally cheaper quarks of 
the Wilson type are used instead, the operator renormalisation is complicated, mainly 
due to the presence of power divergences |18|19j . However, it should be noted that 
power divergences seem unavoidable, even with chirally symmetric regularisations, if 
the charm quark is not included as a dynamical degree of freedom in the low-energy 
approximation to the Standard Model. 

In this paper we advocate the traditional approach of computing K ^ tt tran- 
sitions, in a regime where lowest order chiral perturbation theory can be expected to 
work reliably. For chirally symmetric regularisations this requirement seems prohibitive, 
even in the quenched approximation, as the volume needs to be relatively large^. For 
Wilson-type fermions this problem is significantly less severe, but the renormalisation 
of the four-quark operators is very complicated. Even more importantly, in quenched 
simulations one is again forced to use relatively large meson masses, due to the pres- 
ence of unphysical fermion zero modes. It is the purpose of the present paper to show 
that both the zero mode problem and the problem of operator renormalisation can be 
solved by using twisted mass QCD with four flavours of Wilson-type quarks. This could 
provide an affordable way of approaching the chiral regime by "brute force", down to 

^See, however, ref. |20I21| for an interesting alternative approach, where the effective theory is 
matched in a finite volume very close to the chiral limit. 
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meson masses which are roughly twice the physical pion mass. In this regime the chiral 
expansion can be applied with confidence. 

Since the arguments leading to our results are fairly complex, we have strived to 
present them in a detailed and fairly self-contained way. Moreover, in the conclusions 
we have summarised our findings in the form of recipes which can be applied in future 
numerical computations. Here we only wish to stress the underlying general principles 
which are at play behind the simplification of the renormalisation properties of the 
K ^ TT matrix elements in tmQCD. In this lattice regularisation, a twisted mass term 
is introduced in the Wilson fermion action. The ratio between the new (twisted) and 
the standard mass terms, expressed by "twist angles", needs to be constrained in order 
to recover continuum QCD from the tmQCD regularisation. We are free to fix the twist 
angles to any value, but once the choice is made it is binding, in that it effectively labels 
a specific mode in which chirality is broken by the Wilson term. This in turn implies 
that a given operator is in general expressed in the twisted theory by a combination 
of all the components of its chiral multiplct. Since chiral symmetry is broken in the 
bare theory, each of these components has different renormalisation properties. Thus, 
judicious choices of the twist angles allow us to map the original operator into one of 
its chiral partners which has better renormalisation properties. This is exactly what 
happens in the case of the operators characteristic of K ^ ir transitions. We are able 
to map the original operator, which in the standard Wilson fermion regularisation has 
a quadratic divergence, into one which only has linear divergences. Since quadratic 
divergences are known to be prohibitively hard to subtract in non perturbative com- 
putations, while linear ones are well within our control, this is a significant gain. This 
softening of the divergence comes at a modest price: the twisted theory is characterised 
by soft symmetry breaking; in particular parity and isospin (in two directions of isospin 
space) are lost. Although these symmetries are recovered in the continuum limit, their 
breaking in the bare theory allows for more counterterms to come into play. In the case 
under consideration, the breaking of parity implies that one extra linearly divergent 
counterterm needs to be subtracted; i.e. an extra renormalisation condition must be 
imposed. This condition amounts to the restoration of parity as the continuum limit is 
taken. These considerations are valid in general. In more specific cases (choice of twist 
angles, improved action, quenched strange and charm quarks) we have shown that it 
is possible to avoid the linear subtraction altogether, by working with a pion source at 
large-time separations. 

The layout of this paper is as follows. After a brief review of the effective weak 
Hamiltonian of the Standard Model, and the effective chiral theory at low energies 
(sect. 2), we recall the renormalisation properties of the relevant four-quark operators 
both in a chirally invariant regularisation and with standard Wilson quarks. In sect. 3 
we introduce twisted mass QCD (tmQCD) with four quark flavours and establish the 
relations between its correlation functions and those of standard QCD. Renormalisation 
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and practical aspects of lattice regularised tmQCD are discussed in sect. 4. We are then 
prepared to describe strategies of how to extract K — > tt matrix elements (sect. 5) and 
we end with a short summary of our findings (sect. 6). This work has appeared in 
preliminary form in refs. j22j . 

2 The A/ = 1/2 rule and lattice QCD 

In order to define the general framework we shortly review some aspects of the AS = 1 
effective weak Hamiltonian of the Standard Model, and its chiral effective theory in 
terms of kaons and pions. 

2.1 Effective weak Hamiltonian 

At scales well below Mw, but above the charm quark mass, weak, CP-conserving, 
AS = 1 interactions can be described with an effective Hamiltonian of the form^ 



Weff = KdK"s^[C+(M)0+(/i)+C_(^)Oj^(^)]+h.C. 



(2.1) 



Here 0^{p) are the four-fermion operators 



^ c 



(2.2) 



renormalised at a scale fi (the subscript R indicates renormalised operators), C±{n) 
are the corresponding Wilson coefficients, and 7^ = ^^{l — 75). In principle, a further 
contribution to the weak Hamiltonian comes from integrating out the top quark. This 
has been neglected here, as it is suppressed by a factor (VtdV^*) / (VudVus) — O(10~^). 

The operators can be classified according to isospin symmetry as follows: 0~ 
is in the 1=1/2 representation, whereas O"*" is the sum of two contributions. 



0^ 



(9+ -I- (9+ 



(2.3) 



corresponding to the representations with 1 = 1/2 and 7 = 3/2 respectively. These are 
given explicitly by 



O 



1/2 



•^3/2 



1 



(S7^n)(^z7^d) + isjf:d)iujS) + 2{s^^MdlU) 



(^7uc)(c7^(i) + (s7^(i)(c7^c) 



(S7^n)(^z7^d) + (S7^d)(t27^^z) - {s^U){diU) 



(2.4) 
(2.5) 



''We follow the standard convention and include efltective quark bilinear operators as counterterms 
in the renormalised four-quark operators. 
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The AI = 1/2 rule implies that the contribution of the 1 = 1/2 component of the 
effective Hamiltonian to if — > tttt matrix elements is much larger than that of its / = 3 /2 
part. This contradicts naive expectations and results in the 1/Nc expansion (cf. \2'6\). 

Due to the chiral structure of the weak interactions, the operators or the effective 
Hamiltonian transform as singlets under the SU(3)r part of the chiral flavour group. 
Their classification according to the irreducible octet and 27-plet representations of 
SU(3)l can be found e.g. in ref. [20]. Here we note that we will be interested in 
regularisations with Wilson quarks and non-standard mass terms, which partially break 
the chiral flavour symmetries. The decomposition of the renormalised operators into 
irreducible chiral representations is then non-trivial and hence of limited use. 



2.2 Effective chiral theory of kaons and pions 

At low energies kaons and pions are the relevant degrees of freedom, and their inter- 
actions can be described by chiral perturbation theory (xPT) j24j . in terms of the 
(Euclidean) action 

d^x— {Tr (d^Ud,,U^^ - Tr (Ux^ + xU^) } • (2.6) 

Here, F is the the pseudoscalar meson decay constant, normalised so that i<V = 
92.4 MeV is the experimental value determined from pion decay. The SU(3) matrix 
field U{x) collects the pion, kaon and eta fields, parameterised as 



U{x) =exp{V2i^{x)/F), (2.7) 

with 

$=( IT- -7r7V2 + r//V6 |. (2.^ 



7r^/V2 + rj/V6 


7r+ 


K+ 




-7r7V2 + r//y6 




K- 


^0 


-27]/y/6 



A chiral flavour transformation {gLidR) £ SU(3)l x SU(3)^ transforms the [/-field 
according to 

U - QlUq^, (2.9) 

and the action is invariant under this symmetry, provided the spurion field x transforms 
in the same way as U . The symmetry is broken explicitly by setting 

X = 2BM, M = d\a.g{'mu,md,ms), (2.10) 

where i? is a constant. Assuming isospin symmetry, m„ = rrid = rh, and expanding the 
[/-field in powers of <I>, one finds the lowest order relation between meson masses and 
quark masses, 

rn^ = 2Bm, rnj^ = B{nis + ifi), (2-11) 
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which imphes, 

X = diag{ml,ml,2m'j^-ml). (2.12) 

The weak vertices for AS = 1 processes which correspond to the effective 4-quark op- 
erators in the Standard Model can now be identified by their transformation behaviour 
under the chiral flavour symmetry. Following the conventions of |2nj we have 



= 2\/2GFKdK*s |^527027(x) + 2gs08{x) + 2g'Mx)^ + h.c. . (2.13) 



Using the notation 

2 

the operators are explicitly given by 



L^, = ^d^UU\ (2.14) 



021 = ^(iM)23(^M)ll + ^(^M)2l(^M)l3, (2.15) 

1 ^ 

= i^Y.^L^)2k{L^)kz, (2.16) 
k=l 

O', = ^F\Ux^ + xU^)23. (2.17) 

Besides the parameters F and B, first order for electroweak processes of AS" = 1 
contains the undetermined couphngs §27, gs and g'^. Moreover, it can be shown that 
the so-called weak mass term proportional to g'g does not contribute to physical kaon 
decay amplitudes |lUI25j . 

By matching experimental results for K tttt decays to the rates obtained at 
lowest order chiral perturbation theory one may obtain a phenomenological estimate of 
gs and (727 • Taking into account experimental information about the scattering phases 



between the pions in the final state, the authors of refs. j23l26j obtained the values 

b8|«5.1, |(727|~0.29. (2.18) 

The clear hierarchy between the couplings reflects the A/ = 1/2 rule in the framework 
of leading order chiral perturbation theory. 



2.3 K ^ IT amplitudes 

An explanation of the AI = 1/2 rule in this framework requires the determination 
of the couplings gs and g'27 directly from QCD. Instead of a direct determination of 
the physical K -kit transitions, it has become customary to study the unphysical 
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K ^ IT amplitudes JOl- Using standard conventions (see e.g the matrix element 
for 7r+ transitions in the effective theory is 

(7r+,q|W^(0)|i^+,p> = V2GFVudV:,F^ | ^527+58) p-q + 2mlg's^ , (2-19) 

where 

p-q = EK{p)EM-p-^, Exip) = ^m^ + p2. (2.20) 

In order to determine the combination |g27 + 58 we note that this term is proportional 
to the product of spatial momenta p • q, while the unwanted g'g term is part of the mj^ 
contribution. One could thus isolate the former term by e.g. evaluating the 
matrix element in lattice QCD for two different pion energies (i.e. two different sets of 
spatial momenta). A technical disadvantage of this approach (as opposed to the one 
adopted with domain wall fermions in [3]) is that it requires computation of lattice 
correlation functions with non-zero spatial momenta, which are noisy. 

In order to obtain both g27 and gs independently one then needs a second matrix 
element. If isospin symmetry is unbroken on the lattice, one may use the decomposition 
(|2.3|) and the corresponding one in the effective theory, where the octet operator only 
mediates transitions with AI = 1/2, whereas the 27-plet operator contains both, AI = 
1/2 and 3/2. Decomposing this operator 



with 



0l/2 
•-^27 



027 = Oli^ + Oli\ (2.21) 



^ [{L^)2i{L^)i3 + (iM)23 {4(^.^)11 + 5(^^)22}] , (2.22) 
^ ^(L,,)2i(^m)i3 + (^m)23 {(^m)ii " i^M] > (2-23) 
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it is easy to check that the — > vr"^ matrix elements of these operators yield two 
independent combinations of gs and g27- If isospin is not a lattice symmetry the isospin 
decomposition of the lattice operators poses a difficult renormalisation problem by 
itself. Rather than addressing this problem, one may instead look at the tt^ 
matrix element, 

(7rO,q|W^(0)|i^°,p> = GFKdK*si"' {(527 -g8)p-q- 2mlg's} , (2-24) 

which yields the combination g27 — gs- However, we anticipate that this matrix element 
may be difficult to evaluate in twisted mass lattice QCD (cf. sect. 5.6). 

Finally we recall the fact that (727 can be related, via chiral symmetry, to Bk in the 
chiral limit. Using an independent determination of Bk as input, one may subsequently 
obtain g^ by studying the vr"'" matrix element of eq. 1)2. 19(1 . 
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2.4 Renormalisation of the four-quark operators on the lattice 

In practice the renormaUsation of the four-quark operators has been one of the major 
stumbhng blocks for lattice regularisations which do not preserve chiral symmetry, 
such as lattice QCD with Wilson type quarks |2Sj. Decomposing the operators into 
parity-even and parity-odd parts, 

= ^vv+AA ~ ^VA+AV (2.25) 

we see that parity ensures that K — > vrvr matrix elements receive contributions only 
from OyA+AV '^l^il^ i^T — > vr matrix elements arise only from Oyy_(_^^. 

To illustrate the importance of chiral symmetry we briefly recall the counterterm 
structure in the case of an ideal regularisation with cutoff 1/a, in the sense that the 
continuum symmetries are preserved (see e.g. |29|). In this case one finds that the 
operators renormalise as follows: 

(Or)va+av = C>VA+AV + ("^c - ^l) i^rus - nid) s-f^d , (2.26) 



/VA+AV 



(Or)vv+aa = ^ [Oyy+aa + - O i^s + ma) sdj . (2.27) 

Here, the multiplicative renormalisation constants and the additive renormalisa- 
tion constants are the same for both operators. Counterterms which vanish by the 
equations of motion have been omitted. This is justified as long as only on-shell corre- 
lation functions are considered. Note that lattice regularisations with Ginsparg- Wilson 
fermions provide explicit examples for the ideal regularisation, if one disregards the 
breaking of the continuous space-time symmetries. Therefore, with Ginsparg- Wilson 
fermions, the renormalisation pattern of the operators can indeed be cast in the form 
(|2. 2612. 27)1 |16j . Furthermore, one expects that the operators are on-shell 0(a) im- 
proved, i.e. leading cutoff effects in their matrix elements are of O(a^). 

In contrast, for Wilson-type quarks one finds the counterterm structure (see e.g. jl91 

(Or)^a+av = ^^[o^A+AV + 4^75^^], (2.28) 

5 

(Or)^v+aa = Z^[j2z^O^ + 4 sd + dtsa^,F^,d\. (2.29) 



k=l 



Here k G {VV -|- AA, VV — AA, SS, PP, TT} labels the parity even Dirac structures 
of the four-quark operators in standard notation [29j . and -F^,y denotes some lattice 
version of the gluon field strength tensor. Once more, counterterms which vanish by 
the equations of motion have been neglected. 
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A first observation is that the parity even and odd operators are renormalised 
differently. Compared to the chirally symmetric regularisation one also finds that some 
of the additive renormalisation constants are divergent as a — > 0: 



This is due to the loss of a factor (mc+m-u) and, for the parity even operator, of a further 
factor (nis + rrid)- Instead of a finite counterterm one now has linear and quadratic 
divergences for the parity odd and even operators respectively. In the latter case there 
is also a subleading finite counterterm. An additional complication arises in the parity 
even operator due to the mixing with four other operators of dimension 6. Finally, 
we note that the renormalised operators are not 0(a) improved. These complicated 
renormalisation properties of the parity even operator have prevented the computation 
of K ^ TT transitions with Wilson fermions |28j . On the other hand, the situation 
for the parity odd operator is not much worse than the case of a chirally symmetric 
regularisation: only a single additive counterterm needs to be subtracted (besides the 
multiplicative renormalisation). This motivates us to investigate whether twisted mass 
QCD may alleviate the lattice renormalisation problem, in a similar way as discussed 
for i^TT and the matrix element for K^-K^ mixing in refs. 130131132] . 

3 Twisted mass QCD with four quark flavours 

Twisted mass QCD has been designed to eliminate exceptional configurations in (par- 
tially) quenched lattice simulations with light Wilson quarks 30 . In its original for- 
mulation, it describes a mass-degenerate isospin doublet -0 of Wilson quarks for which, 
besides the standard mass term, a so-called twisted mass term iiJ,q'ip^5T^tp is introduced. 
The properties of tmQCD have been studied in detail in [SJ, where, in particular, its 
equivalence to standard two-fiavour QCD has been established. We discuss here the 
extension of this framework to four quark fiavours. 

3.1 Classical action 

We start by considering tmQCD in the continuum with the fermionic action 
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- {rric - niu) (m^ - nid) 




(2.30) 




(3.1) 
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where, ifr = {u,d,s,c). The standard and twisted mass matrices have the form, 
respectively: 

m = diag (Mu cos a, Mrf cos a, Ms cos /3, Mc cos /3) 

= dia.g{mu,md,ms,mc) , (3.2) 

/I = diag (M„ sin a, —M^; sin a, Ms sin /3, —Mc sin /3) 

= diag(^„,/id,/is,/ic) • (3.3) 

Hence the theory has six independent parameters, namely the four radial quark masses 
Mi {i = u, d, s, c), with = mf + and the two twist angles a, /3. In other words, 
the four standard mass parameters rrii and the four twisted mass parameters are 
constrained by 

tana = = , tanp = — = . (3.4) 

m„ rud nis rUc 

This framework is a natural extension of two-flavour tmQCD, with the property that 

the quark mass terms remain flavour diagonal. At vanishing twist angles a and /3 one 

recovers the standard QCD action of four quark flavours, while for /3 = and M„ = 

Md the two-flavour version of tmQCD of ref. (Slj is reproduced, with two additional 

untwisted quark flavours s and c. 

The action (|3.1j) is form invariant under the axial field transformations 

^ ^ ^' = R{a, /3)V', ip^ip' = ^R{a, (3), (3.5) 

with 

R{aJ) = exp I (arf + /3r^) I , (3.6) 

rf = diag(l, -1,0,0), (3.7) 

ri = diag(0, 0,1,-1). (3.8) 

The new mass term is given by 

m' + f75/x' = R{2a, 2(3) [m + ^75/1] , (3.9) 

and can again be parameterised as in eqs. (|3.2I3.3|) . The corresponding twist angles 
{a', 13') = {a — a, (3 — P) are related to the mass parameters m',/x' as in eq. ()3.4|) . 
The radial masses Mj remain invariant. With the choice a = a, P = P one obtains 
n' = and m' = diag(Mu, M^, M^, Mc), i.e. four-flavour QCD with a standard mass 
term. This demonstrates the equivalence of tmQCD and standard QCD at the level of 
the classical action. It also implies that both theories share all the symmetries. For 
example, a parity transformation, x ^ x = (xq, — x), is realised in tmQCD by the field 
transformations 

V'(x) ^7oi?(2a,2/3)V(x), V^(x) ^ ^(x)i? (2a, 2/3) 70, (3.10) 
with angles a,/3 given by eq. ()3.4|) . 
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3.2 Relations between composite fields 

In general, the axial transformation (|3.5() induces a mapping between composite fields. 
For quark bilinear operators, it can be written in a compact form: 

ij.Tij, = cos((/)r)V5^rV';+ir/rsin((^r)V5^r75V'-. (3.11) 

Here we have assumed that the spin matrices T either commute or anticommute with 
75. The phases are then defined by 

,r / +1, if{75,r} = 0, 

and 

w« = -Ud = ^s = -Wc = -. (3.13) 

These relations establish a "dictionary" between composite fields in theories with dif- 
ferent parameterisations of the quark mass terms. As shown in ref. [HJI the very same 
relations hold between properly renormalised fields in the quantum theories, provided 
the renormalisation scheme respects the chiral and fiavour symmetries of the massless 
continuum theory. 

For illustration and later use we quote a few specific examples. With the notation 
for quark bilinear fields 

Sij = ipiipj, Pij = i^a^i^j, A^^ij = ipa^jaTpj, v^^ij = ipa^,ipj, (3.i4) 

we find, for instance 

Sus = cos(^)5L-isin(^)p^„ (3.15) 

Pus = cos(^)p^,-zsin(^)5;„ (3.16) 

A^,us = COS A'^^^, + ism [^) V^^^„ (3.17) 

V^,,us = cos (^^) V^,^^, + isin (^^) A'^^^^. (3.18) 

We observe that the axial transformation amounts to a rotation by the angle (a + /3)/2 
in the coordinates (Pus, —iSus) and by the angle (a— /9)/2 in the coordinates {Aus,iVus)- 
With a sign flip a —a the same equations hold if the up quark is replaced by a down 
quark. As for operators with up and down quarks only, we find 

Pan = PL Al = A'l, (3.19) 
where we have used the isospin notation 

4 = ^ (\,uu - A^,,dd) ■ (3.20) 
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3.3 Partial current conservation 

With non-degenerate quark masses all chiral and flavour symmetries are broken explic- 
itly. This is expressed by the non-vanishing r.h.s. of the PCAC and PCVC relations. 
For the flavour non-diagonal operators with i 7^ j, they read, 

diiA^^ij = (rrii + rnj)Pij + + fij)Sij, (3.21) 
d^iV^,ij = (mi - mj)Sij + - fJ.j)Pij. (3.22) 

These relations take their usual form when expressed in terms of the axially transformed 
fields. Indeed, setting a = a and j3 = [5, the PCAC and PCVC relations become 

a^^;,,, = (M, + M,)/^^, dX,j = {M, - M,)S[^. (3.23) 

3.4 Mapping of four-quark operators 

Using eq. (|3.1H) we obtain the relations for the four-quark operators of interest: 

Ovv+AA = cos(^)o;±+AA + isin(^)o;±_,AV> (3-24) 

O^A+AV = cos(^)o;±+Av + ^si^(¥)oVWA- (3-25) 

The transformation is a rotation by the angle (d -|- f3)/2 in the plane spanned by 
(Oyy_i_^^, iOyj^j^^) aud can be easily inverted by reversing the sign of the twist angles. 
Setting again a = a and /? = /? it is worth noting two special cases: 

(i) a = /3 = 7r/2: with this choice the operators are directly interchanged, as 
eqs. (|3.24I3.25|) become 

^VV+AA = ^^VA+AV ^VA+AV = ^^VV+AA" (3.26) 

In terms of the quark masses these twist angles correspond to vanishing standard 
mass parameters, 

ruu = rud = rus = rric = 0, (3.27) 
and positive twisted mass parameters for up and strange quarks. 

(ii) a = —P = Tr/2: in this case one finds 

^VV+AA = ^VV+AA' ^VA+AV ^ ^VA+AV (3.28) 

In terms of quark masses the only difference with respect to the previous case is a 
change of sign for the twisted mass parameters of the strange and charm quarks. 

We will refer to both these cases as "fully twisted", since the physical quark masses 
Mi are then determined by the twisted mass parameters alone. In the fully twisted 
cases the operator OyA+AV ^™ P^^y either the parity even or the parity odd 

operator, depending on the sign of /?. 
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3.5 Relations between correlation functions 



Consider, in the formal continuum theory, Euclidean correlation functions of the form 

{0[^,i^]) = Z-^ [ 0[V',V^]e-^ (3.29) 

where 0['i/','i/'] denotes some multilocal gauge invariant field. The fermionic part of the 
action is given in (|3.H) . Performing the field transformation eq. (|3.5j) and changing 
integration variables in the functional integral one derives the identity 

(O [R{aJ)^,i}RiaJ)])^^^^ = {O [V', V^] >(,_5,^_^) (3-30) 

where the integration variables are ■0 and ip on both sides of the equation, and the 
Jacobian is unity, due to 

det fi?(a,/3)) = 1. (3.31) 



Setting a = a and /? = /3 we thus obtain identities between standard QCD and twisted 
mass QCD correlation functions, which we can use in two alternative versions: 

(0[V,V5]>(„,^) = (0[i?(-a,-/3)V,V^i?(-a,-/3)])(o_o), (3.32) 
<0[V,V5])(o,o) = (0[i?(a,/3)V',V^i?(a,/3)]>(„^^). (3.33) 

Hence, a given correlation function in tmQCD with twist angles (q, (3) is interpreted 
as the linear combination on the r.h.s. of eq. (|3.32j) . If instead we are given a standard 
QCD correlation function, then eq. (|3.33|) tells us how it is represented in tmQCD at 
twist angles (a,/?). The dictionary just established for composite fields can be used 
for the integrands in eq. (|3.32j) . provided ijj' and ip' are identified with the integration 
variables on the r.h.s.. As shown in ref. the relations derived for the formal con- 
tinuum theory are realised between the renormalised quantum field theories, provided 
the renormalisation procedure is set up with some care. 



4 Lattice regularised tmQCD 

To regularise tmQCD on the lattice, we assume that the gauge fields are represented 
by some standard lattice action (e.g. Wilson's plaquette action), whereas quarks are 
described by the lattice action 

5f = a^^^(x)(L'iy + m + i75/x)'0(x), (4.1) 

X 

with D\Y being the standard, possibly 0(a) improved Wilson-Dirac operator (for un- 
explained notation see |S3])- In order to ensure the equivalence between tmQCD and 
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standard QCD, the bare mass parameters, rrii and //j (i = u, d, s, c), will be subject to 
constraints. However, these are now only implicitly defined by eqs. (jSHJ, as they will 
be imposed on the properly renormalised mass parameters. Before discussing renormal- 
isation of lattice tmQCD we consider the possible choices for the bare mass parameters 
leading to feasible numerical simulations. 



4.1 Quark determinant and practical choices of parameters 

Whether a given lattice action is suitable for numerical simulations obviously depends 
on the available algorithms. The answer may therefore change with time. At present, 
the starting point for most algorithms is the functional integral after integration over 
the quark fields. The resulting determinant is then part of the effective gauge field 
measure and is generally required to be real and positive. In the case of lattice tmQCD 
with four quark flavours as introduced here, each flavour contributes a factor to the 
determinant. Assuming the light doublet to be mass degenerate, 

mu = md = mi, fj,^ = -fj^ = fj^^ (4.2) 

the determinants of the two light flavours can be combined to yield 

det [Dw + + ifin^T^) = det {Dw + m^)^ {Dw + m) + /^F , (4.3) 

where the determinant on the r.h.s. is taken in a single flavour space .31 . Hence 
the determinant of the light doublet is positive at non-zero fj-i, irrespective of the 
background gauge field. Integrating over strange and charm quarks one obtains 



det 



{Dw + rus)^ {Dw + rUc 



+ i ^icTo {Dw + ms} + i iisl5 {Dw + "id , (4.4) 

which, in general, is not a real quantity. Discarding the unphysical situation of mass- 
degenerate strange and charm quarks, the only way to ensure the reality of the deter- 
minant is to employ untwisted strange and charm quarks, /i^ = = 0, as the fermion 
determinants for individual Wilson quark flavours are real. However, positivity of the 
determinant is then not guaranteed, unless the quark masses are large enough. Detailed 
statements depend on the lattice gauge action and the chosen simulation parameters, 
but for practical purposes the strange quark mass might be heavy enough (see, e.g. |3l]). 

In view of this situation we distinguish the following practical options for lattice 
tmQCD simulations: 

• Quenched simulations or partially quenched simulations with A'^f = 2 dynamical 
light and mass degenerate quarks. Strange and charm quarks remain quenched 
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and do not contribute to the determinant, which is therefore guaranteed to be real 
and positive (for 7^ 0). Hence, one has almost complete freedom in the choice 
of the mass parameters and therefore of a and f3. In the quenched approximation 
one would just require a non-zero twist angle a in order to ensure the absence of 
unphysical zero modes. 

• A'^f = 3, 4 dynamical quark flavours with two mass degenerate light quarks, and 
untwisted strange and charm quarks. This implies the choice for the second twist 
angle /? = 0, while a remains unrestricted. 

Note that these considerations pertain to our specific choice of flavour diagonal mass 
terms. If one allows for flavour non-diagonal mass terms one may indeed be able to 
obtain positive and real fermion determinants for any combination of non-zero quark 
masses, including non-degenerate light quarks [SHI. However, this renders the relation 
between standard QCD and tmQCD correlation functions more complicated, and the 
flavour structure has to be dealt with explicitly in the numerical computation of quark 
propagators. 

4.2 Renormalisation of lattice tmQCD 

Renormalisability is based on general properties of Quantum Field Theory, such as 
locality and unitarity, which may not hold in the (partially) quenched approxima- 
tion of QCD. In this paper we do not try to solve this general problem. Instead we 
take the naive point of view that the (partially) quenched theory is well-defined and 
can be obtained by suppressing the fermion determinant with respect to some or all 
flavours. This is assumed to incur no drastic change to the theory's renormalisability. 
While we discuss the counterterm structure for the full theory with A'^f = 4 dynam- 
ical quarks, we occasionally indicate modifications which should be expected in the 
(partially) quenched case. 

4-2.1 Symmetries of the lattice action 

When the mass terms are taken to vanish, the tmQCD lattice action reduces to the 
usual (possibly 0(o) improved) Wilson action for A'f = 4 quarks. Hence, the massless 
tmQCD lattice action (m = /x = 0) is invariant under U(4) flavour rotations. 



(4.5) 



as well as under discrete transformations, such as parity (x — > x = (xq, — x)) 



V'(x) 7o'0(*) 



Tp V'(^)7o 



(4.6) 
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and charge conjugation 

V'(x) ^ C-^i^ixf, i) -ipixfC (4.7) 

(the charge conjugation matrix satisfies 7* = —Cj^C~^).^ 

Once the mass and twisted mass terms are added, some of these symmetries are 
lost. The introduction of the m term imphes loss of the f7(4) flavour symmetry, unless 
m oc 1. Upon introducing the fx term this symmetry is also lost. A reduced symmetry 
survives however, if /x cx diag(l, — 1, 1, — 1). This consists of isospin rotations of the 
SU{2) subgroups, along the generators rf and defined in eqs. (|l-i.7|) and (|;-i.<S|) . Charge 
conjugation is unaffected by the mass terms, while parity is clearly broken by the twisted 
mass term.^ 

Upon adding the two mass terms, one may determine the counterterm structure 
by treating m and fi as spurion fields which transform under these symmetries so that 
the massive theory remains formally invariant. In particular, under a U(4) flavour 
transformation (|4.5j) one assumes that 

m^VmV\ fi^VfiV\ (4.8) 

while invariance under the parity transformation (|4.6|) is obtained by assuming 

m — > m , A* ^ — A* ■ (4-9) 

However, we stress again that in tmQCD with Wilson quarks and nonvanishing twisted 
masses physical parity, as defined in eq. (|3.1fl|) . is broken at 0(a), and is recovered only 
in the continuum limit. Analogous considerations hold in the case of the U(4) vector 
flavour symmetry. 

It is sometimes useful to restrict attention to a subset of the general flavour trans- 
formations (|4.5|) . such as flavour exchanges 

ipi^ipj, ipi^i'j, (4.10) 

for which the spurion transformations (|4.8|) reduce to an exchange of the corresponding 
mass parameters, 

rrii ^ nij, f^j. (4-11) 

®We have not indicated the parity and charge conjugation transformations of the gauge field, as 
they are never used in this work. 

^The reader should not be misled by the fact that the redefined parity transformation 13.1011 is a 
symmetry of the classical action. Although the fj, term is invariant under this symmetry, the Wilson 
term is not. 
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4-2.2 Renormalised parameters 

The physical interpretation of tmQCD hinges on the knowledge of the twist angles in 
the renormalised theory. As these are determined by ratios of the renormalised standard 
and twisted mass parameters we first discuss renormahsation of the bare parameters in 
the action. 

The counterterms to the action of dimension < 4 are, 

triF^^F^^}, ipTp, tr(m)V;V', tpmip, ztr(/x)V'75''A, ^^/^TsV", (4.12) 

where Ff^i, denotes the gluon field strength tensor. The first counterterm implies a 
multiplicative renormahsation of the bare gauge coupling. As for the quark mass renor- 
mahsation the situation is complicated by the trace terms, which would be absent in a 
chirally symmetric regularisation. For standard Wilson fermions the situation has been 
discussed in refs. |,S6p,S7j . which can be summarised and extended as follows: first one 
decomposes the mass matrices in non-singlet and singlet pieces, 

15 15 

m = ^m''A'' + m°l4, /x = ^ //°A" + (4.13) 

a=l a=l 

where A"^, a = 1,...15 are the generators of SU(4). This decomposition applies to 
general Hermitean 4x4 matrices, but with our choice of flavour diagonal mass matrices 
only 3 of the non-singlet coefficients are non-zero. These correspond to the 3 diagonal 
and traceless generators of SU(4) and are therefore linear combinations of (twisted) 
mass differences. The singlet coefficients are given by 

mO = ^tr(m), /x^ = itr(/x) , (4.14) 

The non-singlet terms are renormalised by the inverse renormahsation constant of the 
non-singlet scalar and pseudoscalar densities, Sij and Pij {i / j) respectively, 

= Z^'m\ /x^ = Zp V^. (4.15) 

The renormalised singlet masses m? and fi^ are then given by 

ml = Z-,\m^ - mcr), /i^ = ZpoV°, (4.16) 

where Zgo po are the multiplicative renormahsation constants of the singlet densities 
= Yl^=i ^ii = Yl!i=i ^ii respectively, and mcr denotes the "critical mass" . 

From this renormahsation pattern we obtain for quark mass parameters of indi- 
vidual flavours 

mR,i = [mi - mcr + (rm - 1) ("2-° - rucr)] , (4.17) 

^^R,^ = Zp 1 [/ii + (r^ - 1)/] , (4.18) 
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with 



ZsZgV, = Z^Z~l (4.19) 



Here, the ratios r^, are finite functions of the gauge couphng, which can in principle 
be determined by axial Ward identities in the chiral limit. A few comments are in 
order, in relation to the practical considerations made in subsect. 14.11 

• The above counterterm structure applies to the complete theory with A'^f = 4 
dynamical quarks, for which our preferred choice of twist angles is q = 7r/2 and 
/3 = (partial twist). The latter condition translates to renormalised twisted 
parameters 

MR,s = m,c = , (4.20) 

which is achieved by setting 

= = . (4.21) 

Due to Hu = — = this implies = Q and thus no flavour mixing in the 
twisted mass renormalisation pattern expressed by eq. ()4.18|) . On the other hand, 
the condition concerning the twist angle a translates into 

mR,/ = , (4.22) 

which, due to eq. (|4.17|) implies that the bare mass of the light quarks must be 
tuned as follows, 

{I — j-\ 

mi = nicr + , . {ms + rUc - 2mcr) , (4.23) 
2(1 + rm) 

i.e. there is an offset of 0(1) to the linear divergence cancelled by rricr- 

• The ratios and depend implicitly on the number of dynamical quark 
flavours. In particular, their quenched values are = = 1, i.e. the flavour 
mixing terms in eqs. (|4.17|) and (|4.18j) are absent in this case !^|. Hence, the 
fully twisted cases a = ±/3 = 7r/2 are obtained by just setting all standard bare 
masses equal to the critical mass rricr- 

• In the partially quenched case it is obvious that valence quark masses do not 
participate in the renormalisation of other quark flavours. More specifically, for 
the case of interest with Nf = 2 light dynamical flavours, the mass renormalisation 
pattern reduces to 

rriR^l = Z^^Tm {mi - m^r) , (4.24) 
m-R.s.c = [ms,c - rricr + (rm - 1) {rni - mcr)] , (4.25) 

= V* ■ (4.26) 
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This implies that the fuhy twisted case, a = ±/3 = 7r/2, requires a quenched-hke 
tuning mi = rus = rric = nicr, while for the partially twisted case, a = 7r/2, f3 = 0, 
the requirement is m; = rricr and fj,s = iJ^c = 0. 

4-2.3 Determination of the twist angles a and (5 

In practice, the twist angles are defined by ratios of the renormalised quark masses 
appearing in the renormalised PCAC and PCVC relations |nj|. In general, one therefore 
needs to discuss the renormalisation of the currents and densities which appear in the 
PCAC and PCVC relations H3.21I3.22() . If we restrict attention to flavour off-diagonal 
quark bilinear operators (i 7^ j), the symmetries imply a multiplicative renormalisation 
in all cases, 

(^R)M,ij = Zf^A^^ij, (Vr)^,,^ = ZyV^^ij, {PR)ij = ZpPij, {S-R)ij = ZsSij , (4.27) 

where the bare fields are defined as in the classical continuum theory, eq. ()3.14|) . It is 
well-known that the scale independent renormalisation constants and Zy as well 
as the ratio Zp/Zo, can be determined^ by imposing axial Ward identities as normali- 
sation conditions 18 . To obtain the scale-dependent renormalisation constants Zp or 
Zs, it is then sufficient to impose a further quark mass independent renormalisation 
condition on one of the densities. The choice is largely arbitrary and does not affect 
the determination of the twist angles as the renormalisation constant cancels in quark 
mass ratios. Given the renormalised quark bilinear fields, and a choice of bare param- 
eters, the 6 renormalised mass parameters (2 for the mass degenerate light quarks and 
4 for strange and charm quarks) can be obtained by solving a system of 6 independent 
equations, which follow from the renormalised PCAC and PCVC relations for different 
flavour combinations. 

Fortunately the situation is much simpler in the cases of practical interest, namely 
the choices a = ib/3 = tt/2 and {a, (3) = (7r/2,0), which we now discuss in more detail: 

• To obtain a = ib/3 = tt/2 all renormalised standard quark masses should vanish, 
'^R,; = ''TT'R,s = 'm'R,c = 0, which implies that all bare standard quark masses are 
equal. The corresponding critical bare mass parameter m-cr can be obtained by 
requiring 

d^,A^,ud = 0, (4.28) 

for some correlation function of the axial current. Note that the renormalisation 
constant Za is not needed here, and the angles a = tt/2 and (5 = ib7r/2 are 
now fixed up to a sign, which is determined by the sign of the bare twisted mass 
parameters. Moreover, the fully twisted cases are special because the values of 

''The determination of Zv could be avoided by using the point spht vector current which follows 
from the exact flavour symmetry of mass degenerate standard Wilson quarks. 
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the renormalised twisted masses are not needed for the definition of the twist 
angles. 

• The choice /? = is easily realised by just setting fig = Hc = 0, as the counterterm 
structure then implies //r,^s = fi^^c = 0. Setting the twist angle a = tt/2 is then 
still achieved by requiring (|4.28j) . However, the counterterm structure (|4.23jl 
implies that this determination of mi must be repeated every time the parameters 
nis and rric are changed. 

Finally we remark that the definition of renormalised quark masses through the PC AC 
and PCVC relations is convenient, as it can be applied both in the (partially) quenched 
and in the full theory. In this way one may bypass the discussion of the quark mass 
counterterms in the partially quenched case. 

4.3 Renormalisation pattern of OyA+AV 

In view of the AI = 1/2 amplitudes, we are interested in the counterterm structure of 
the operator OyA+AV ™ presence of twisted mass terms. We apply again a spurion 
analysis. In particular, the flavour structure of the operator suggests to consider the 
two flavour exchange symmetries 

d^s, md^nis, fid^fis, (4.29) 
u ^ c, rriu ^ rric, fJ-u ^ fJ-c, (4.30) 

combined with charge conjugation and spurionic parity discussed in sect. 14.2.11 Power 
counting and the behaviour under these symmetry transformations imply the following 
renormalisation pattern: 

(Or)±a+AV = {OvA+AV + 4 Psd + 4 Ssd} ■ (4.31) 

Note that, compared to standard Wilson quarks, the only change consists in the addi- 
tion of the scalar density. The symmetries also determine the behaviour of the coeffi- 
cients in the small a expansion. For the leading terms we find: 

4 = ^\^{mc-mu){ms-rnd)4,a + il^c- ^J'u){^J's- fJ'd)4,b} +^(^)^ (4.32) 
4 = ^ |("ic - "l„)(/is - /id)c|Q + (/ic - - "id)c5^6} + 0(1), (4.33) 

where the coefficients on the r.h.s. are functions of the bare coupling only. We observe 
that both coefficients Cp^ are linearly divergent in general. Simplifications occur in 
special cases: if down and strange quarks are taken to be mass degenerate, = m^. 
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= A*d; then both coefficients vanish identically due to the exact symmetry under a 
charge conjugation combined with the flavour exchange symmetry, s <-> d. Another 
important simplification occurs in the fully twisted cases discussed above, since 

rriu = rud = ms = rric ^ c| = 0(1), (4.34) 

so that one is left with the linear divergence in Cp, which takes the form 

Cp = ~(/^c - ^J'u){^J's - tJ-d)c$b + 0(a), (4.35) 
i.e. subleading terms of 0(1) are found to be absent. 

4.4 0(a) improved action for fully twisted quarks 

While systematic 0(a) improvement of matrix elements of four-quark operators is be- 
yond the scope of this paper, the discussion of the fully twisted case in section 5 shows 
that some simplifications occur if the action is 0(a) improved. To investigate the 
structure of the 0(a) counterterms in the special case of the fully twisted tmQCD ac- 
tion, we make use of results of ref. [SHI) which were obtained for two-fiavour tmQCD 
with degenerate quarks. These results are directly relevant here, since, as discussed in 
subsection 4.1, a full twist of all quarks is implemented either in the quenched approx- 
imation or in the partially quenched case with dynamical (and mass degenerate) up 
and down quarks. 

In the following it is assumed that the massless theory has been 0(a) improved by 
the addition of the Sheikholeslami-Wohlert term to the action and of the counterterm 
proportional to ca to the axial current [HS]- Let us now suppose that the fully twisted 
case has been realised by tuning the standard bare masses to their critical value as ex- 
plained in subsection 14 . 2 . 31 The 0(a) effects in the renormalised light mass parameters 
can then be inferred from the result of ref. |HSI, as the quenched strange and charm 
quarks may not generate any additional counterterms: 

"iR,/ = Z^obrmafif + O(a^), (4.36) 
I^K^i = ZpVi + 0(a2). (4.37) 

To find the counterterm structure for the strange and charm mass parameters, we use 
the fact that both quarks remain in the quenched approximation. Hence, one may 
imagine that either quark is a member of a mass degenerate quark doublet, with part- 
ners that do not enter any of the correlation functions, and are therefore completely 
decoupled from the theory. The renormalisation for each doublet then proceeds inde- 
pendently with the result (i = s,c), 

m,i = ^pVi! (4-39) 
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where the new coefficients 6^, are non-zero only if the hght quarks are dynamical. At 
this point we recall that in the two flavour theory the set of 0(a) mass counterterms 
as determined by the lattice symmetries is over-complete j38j . More precisely, an axial 
rotation re-parameterises the mass counterterm basis, implying that there is a one- 
parameter family of equivalent 0(a) improved theories. This freedom may then be 
used to set one of the improvement coefficients to a fixed value. Hence, in the fully 
quenched case we may now choose to set 

bmi = bm, = bm, = 0, (4.40) 

and, as a result, all the renormalised quark mass parameters are 0(a) improved. In 
the partially quenched case we may set bmi = 0, and for i = s, c, 

6^, = -b'm^f^f/f^l (4.41) 

We conclude that in the fully twisted cases, there is a choice for the 0(a) mass coun- 
terterms such that the action is 0(a) improved, provided that the action and the axial 
current have been 0(a) improved in the chiral limit. 

5 K ^ ir transitions from twisted mass lattice QCD 

We now show how to extract — > vr"^ matrix elements from tmQCD correlation func- 
tions. We will mostly be concerned with the determination of the additive counterterms 
to the renormalised operator, or with corresponding subtractions of its correlation func- 
tions. The multiplicative operator renormalisation can be performed using standard 
techniques 'Iffl"4() l29p41j and does not pose any further conceptual problems. 

In what follows we will refer to QCD with a standard parameterisation of the quark 
mass term as "standard QCD" , as opposed to QCD with a twisted mass term. Their 
respective renormalised correlation functions are thus related by a chiral rotation of the 
fields, according to the discussion in section |31 

5.1 From correlation functions to matrix elements 

In standard QCD, the 7r+ matrix elements can be obtained from the Euclidean 

3-point function 

^Kni^^ = ('/57r(x)(OR)vv+AA(?/)V5^(^))(0,0) ' (^-^^ 

where (^tt and 99/^ are interpolating fields with the quantum numbers of the charged 
pion and kaon, respectively. Simple examples are the local fields, 

ip-„{x) = Pduix), ^k{z) = Pus{z), (5.2) 
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but more complicated choices are possible and indeed often necessary to obtain a signal 
in numerical simulations. Based on the transfer matrix formalism, the spatial Fourier 
transform 

x,z 

can then be shown to behave, for xq — > oo and zq — oo, as 

G%T,{xo,zo;p,q) ~ Z^Z/^exp(-£'^(p)xo + -E'ii-(q)2:o) 

x(^+,p|(OR)±v+AA(0)|i^+,q>- (5.4) 

Here, neglected terms are exponentially suppressed contributions of higher intermediate 
states in either the pion or the kaon channel. The constants and Zk, as well as 
the pion and kaon energies -E'7r(p) and £'x(q) can be determined from the 2-point 
functions of the interpolating kaon and pion fields, so that one can isolate the desired 
matrix element of the renormalised operator. In fact the matrix element itself is used 
to completely specify the renormalised operator, by imposing the condition 

<7r+,0|(OR)±v+AA(0)|i^+,0> =0. (5.5) 

Note that, even with a chirally invariant regularisation, such a condition may be used in 
order to determine the 0(1) subtraction coefficients in eq. ()2.27|) .^ This arbitrariness 
is closely related to the fact that the i^T — > vr matrix elements are not directly linked to 
a physical observable. Indeed, it is easy to see that the subtracted counterterm does 
not contribute to physical K vrvr amplitudes. This situation is reflected in chiral 
perturbation theory, where the condition 1)5. 5|) determines the unphysical coupling g'^. 
It is actually convenient to impose the renormalisation condition (|5.5|) . as it facilitates 
the matching to chiral perturbation theory. 



5.2 Renormalised correlation function in tmQCD 

In the twisted mass QCD framework, we expect that the properly renormalised corre- 
lation functions are mapped onto their standard QCD counterpart in the same way as 
in the classical continuum theory jS^. According to Sect. 3 we can thus write 

GKnix,y,z) = (^Lpn{x) cos (^^^^ {OB.)y^^j^j^{y) 



-ismi^-^j {OK)yj^^j^y{y)^ipK{z)j^^^y (5.6) 

^Instead of the above condition, in simulations with Ginsparg- Wilson quarks it is customary to 
determine the subtraction coefficient by imposing that the operator matrix element of the unphysical 
K transition (with the Tf-meson at rest) vanishes. 
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where (^^^ and (px denote the transformed source fields, which can be worked out once 
a specific choice of pseudoscalar sources (/?7r and has been made. For instance, to 
reproduce the local source fields of eq. 1)5. 2(1 up to terms of 0(a), one defines 

(^x = cos ( ^ ) (Pr)us + isin ( ^ ) (5r)„s, 'f'-K = {PKldu- (5.7) 



Note that for generic twist angles both parity components of appear in the r.h.s. of 
(|5.6() . while we would like to be left only with the operator O^p^j^p^, due to its nicer 
renormalisation properties. We therefore turn the tables and determine the physical 
interpretation of the tmQCD correlation functions with this operator alone: 

('i?7r(x)(0R)vA+Av(y)^^(^))(^^) = ^99^(x) COS (Or)^A+AV 

isin (^^) {OK)yY+AAiy) 



(0,0) 



= isin(^^f^)G|^(x,y,z). (5.8) 



Here, the last step is trivial in the fully twisted case, a = (3 = 7r/2, for which the 
cosine in the first term on the r.h.s. vanishes. On the other hand, for generic values 
of the twist angles the cosine multiplies a standard QCD correlation function which 
vanishes by parity. Hence, by computing the l.h.s. of eq. (|5.8|) in twisted mass QCD one 
obtains the desired correlation function up to a known factor, and the K ^ tt matrix 
element is then obtained in the standard manner. It is important at this point, however, 
to stress that eq. (|5.8|) holds exactly only in the continuum limit, and assumes that 
physical parity, defined through 1)3. 1U() . is properly restored in renormalised tmQCD. 
In principle this can be achieved in two ways: given the twist angles, it may be possible 
to impose renormalisation conditions so that the renormalised composite fields assume 
the expected continuum behaviour under parity, up to cutoff effects. Alternatively, one 
may be able to identify the parity violating contributions to the correlation function 
on the l.h.s. of eq. 1)5. 81) . so that the parity conserving QCD correlator can be 
obtained. Moreover, this is only required in the asymptotic regime where the matrix 
element is extracted. 

In practice, a combination of both methods seems advisable. Hence we begin by 
identifying three sources of parity violating contributions to the tmQCD correlation 
function in the l.h.s. of eq. (|5.8|) : 

(i) The operator Oy^_,_^y has linear divergences of either parity, cf. eq. (|4.31j) . 

(ii) In general, the interpolating fields (pK,'^-K may break physical parity. However, 
once the twist angles a, (5 are defined from the renormalised PCAC and PCVC 
relations, (pKi'^-w may be chosen so that parity holds up to 0(a) effects (explicit 
examples will be given in the following). 
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(iii) In general, the action will have parity breaking 0(a) counterterms. In the spirit 
of Symanzik's effective continuum theory, these counterterms can be treated as 
operator insertions, and thus contribute at 0(a) to the correlation function. 

As a consequence, unwanted 0(1) contributions which originate from parity violation 
may arise by combining the parity violating linear divergence of the four-quark op- 
erator with parity violating 0(a) terms from either the action or the interpolating 
fields. Hence, if the parity breaking linear divergence of the four-quark operator is 
subtracted by imposing an appropriate renormalisation condition, the parity violating 
contributions to the correlation function are 0(a). 

In the following we make two specific choices for the twist angles a and /? (cf. sub- 
sect. 4.1). For both cases we explain in detail how to obtain the K ^ it matrix elements 
from the tmQCD correlation functions on the l.h.s. of eq. (|5.8p . 

5.3 The fully twisted case a = /3 = tt/2 

For simplicity we specify local source fields which, for a = P = tt/2 read 

(fx = i{S-R)us, ipTT = {PR)du- (5.9) 

We have already pointed out that the fully twisted case is special in that the operators 
of interest are directly matched to each other; thus the K'^ 7r+ correlation function 
(|5.1|) simplifies to 

GK-,ri^,y,z) = {iPR)du{x){OK)^j^^j^y{y){Sjl)usiz))(^/2,7r/2) ■ (5-10) 

It follows that, up to cutoff effects, the renormalised four-quark operators have definite 
physical parity "P, viz. 

7'[(Or)±a+av] = +1 ■ (5-11) 
On the other hand, for the dimension-3 counterterms we have 

V[Psd] = -1 V[S,d] = +1 , (5.12) 

i.e. the linearly divergent coefficient Cp multiplies a parity odd counterterm, while the 
finite counterterm proportional to is parity even. This implies that, out of the two 
renormalisation conditions required to determine the counterterm coefficients Cp and 
c^, the one determining Cp must restore physical parity of the renormalised operator. 
For example, the condition 

((Or)va+av(^) (^R)rf.(y)>(^/2,./2) = 0' (5-13) 
enforces parity conservation, as can be seen from its translation to standard QCD, 

<i(OR)±v+AA(^)(^R)rf.(y)>(o,o) = 0- (5-14) 
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In terms of the bare operators eq. (|5.1.S|) reads 

= (OvA+Av(^) ^'^«(y)>(V2,^/2) 
+ C^{Psd{x) Pdsiy))(^/2,n/2) 
(7r/2,7r/2) • 

(5.15) 

Note that this equation indeed determines the coefficient Cp, as the last term on the 
r.h.s. satisfies 

{Ssd{x) Pdsiy)) (^/2,n/2) = i^pZs)'^ {iSK)sdix) iPR)dsiy))(o,o)' (5-16) 

and thus vanishes by parity up to 0(a). In particular, the prefactor c^, being of 0(1), 
does not affect this conclusion. 

We may now insert the renormalised operator (OR)Y^_(_^y (with Cp determined as 
above and as yet undetermined) in the correlation function (|5.1fl|) . As the mesonic 
source operators of eg . (15.91) are also parity eigenstates, 

V[Sus] = -1 r[Pdu] = -1 . (5.17) 

the remaining parity breaking effects in the correlation function (|5.1flj) are at most of 
0(a), being combinations of 0(a) effects from either the action or the mesonic source 
fields with an 0(1) term of the four-quark operator. It thus remains to determine the 
finite counterterm proportional to c^. This can be done by imposing the condition 
(|5.5|) on the vr"*" matrix element which, in the current tmQCD framework, is 

obtained from the correlation 1)5. 1U|) . 



5.4 Sidestepping the power divergence 



As explained in subsection 4.4, the action in the fully twisted cases can be considered 
0(a) improved, provided 0(a) improvement has been correctly implemented in the 
massless theory. In this situation, the above renormalisation procedure can be further 
simplified, as the determination of the coefficient Cp from the condition 1)5. 15(1 can be 
avoided. 

To see this we write explicitly the correlation in tmQCD: 

G%^{x,y,z) = ZpZsZ^ {Pdu{x)Oyj^^j^y{y)Sus{z))(^/2,n/2) 

+ Cp {Pdu{x)Psd{y)Sus{z))(-^/2,TT/2) 

+ c| {Pdu{x)Ssd{y)Sus{z))(^/2,n/2) (5-18) 



26 



and examine in detail the correlation proportional to Cp on the r.h.s.. The operator 
product PduPsdSus is an odd parity eigenstate, giving vanishing 0(1) contributions. 
The 0(a) corrections (which give rise to 0(1) contributions to the original correlation 
function G^^) may arise from three sources, as listed in subsect. 15.21 

(i) The parity-even 0(a) correction to the operator insertion Pgd has the form a(^s + 
HdjSsd, due to the symmetries of the lattice action. This term amounts to a re- 
definition of and can thus be ignored. 

(ii) The even parity 0(a) corrections to the kaon source Sus are proportional to 
{fJ-u + fJ's)Pus^ while the analogous counterterm for the pion only appears at O(a^) 
(for mass degenerate light quarks); the 0(a) correction to the Fourier transformed 
correlation function is then given by 

a^^e*(P™)(Pd.(x)P,,(0)P„,(z))(^/2,./2) e^o^Ks(q), (5.19) 

x,z 

where we have exhibited the expected asymptotic factor for large negative zq, 
and K'^ denotes the lowest scalar excitations in the kaon channel. It is not clear 
what the exact nature of this state is. Candidates are one-particle states with the 
opposite parity, or multi-particle states. For instance, a J = two-particle state 
consisting of a kaon and a neutral pion shares all the lattice symmetries with the 
kaon, and might therefore be the next lightest state with opposite physical parity. 
In any case, we expect the energies of the excited states to be much higher than 
the kaon energies themselves. Therefore, in the asymptotic regime, the relative 
suppression factor exp{zo [Ej^s{q} — £'x(q)]} of the second term on the rhs of 
eq. (|5.18|) should be significant. 

(iii) Finally, as the 0(a) contribution of the action vanishes with the appropriate 
choice of the counterterm basis (cf. subsect. 4.4), the action contributes at most 
to O(a^). This combines with the linear divergence of the operator to yield an 
0(a) effect to G^^. Here it is crucial that the action be 0(a) improved, as the 
insertion of action counterterms into the correlation function comes with a space- 
time integration, which makes the effect non-local. In particular, it would not be 
possible to argue in terms of the asymptotic regime of the correlation function as 
we did in (ii) with the localised counterterms to the kaon source field. 

We conclude that in the fully twisted case the Cp term in eq. ()5.18|1 need not 
be subtracted explicitly, as it leads at most to 0(1) contributions to which are 
exponentially suppressed, due to the energy gap between scalar and pseudoscalar sec- 
tor. Remarkably, the only subtraction to be carried out explicitly is the finite, parity 
conserving counterterm proportional to Cc. One is thus left with a renormalisation 
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structure that closely resembles the one found in regularisations that preserve chiral 
symmetry. 

5.5 The partially twisted case with /? = 

In this case, the operator Oy^_|_^y is mapped to a linear combination of parity odd 
and parity even pieces (cf. eq. (|,3.25|) ). Therefore, it is not clear how to impose a 
parity restoration condition for the operator. We therefore proceed in two steps: first 
we subtract the two power divergent counterterms with some arbitrary prescription, 
e.g. by imposing the conditions 

((O)^A+Av(^)(^)rf^(y)>(./2,0) = 0' (5.20) 
((0)VA+Av(^)(5)^^(y)>(./2,0) = 0- (5-21) 

Here Oy^_,_^y is the subtracted bare operator, 

Oya+AV ~ ^VA+AV + ^%Psd + c^Ssd (5.22) 

and the coefficients Cp^ are determined by eqs. H5.2UI5.21|) . This bare operator is only 
logarithmically divergent, and its renormalised counterpart is obtained by multiplicative 
rescaling with the usual renormalisation constant, 

(Or)va+av = ^"^^^A+AV (5-23) 

However, such a renormalised operator does not have the physical parity properties that 
ensure the restoration of parity in the continuum limit. Rather, a correctly renormalised 
operator will be of the form 

(Or)va+av = {OvA+Av + ^4Psd + a4s,4 , (5.24) 

with still unknown, but finite coefficients Acp^. These should be determined so that 
the hadronic matrix element of the renormalised operator does not contain unwanted 
contributions of the wrong physical parity. This does not necessarily determine the 
renormalised operator itself. Indeed, it turns out that only a single coefficient needs 
to be determined, the other giving rise to contributions which are either of 0(a) or 
exponentially suppressed or both, depending on the choice of interpolating kaon and 
pion fields. 

To see this more explicitly, we have a closer look at the additive counterterms, Pgd 
and Ssd- It is convenient to split both fields in pieces which are either even or odd 
under the physical parity transformation. 

Pds = {Pds)even ^ (^i^s)odd ' ^ds = {Sds)even + ('^'is)odd • (5.25) 
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From the correspondence between the renormaUsed fields in tniQCD and standard QCD 



(cf. subsection 3.2) one then infers the explicit expressions, 

(^ds)even = sin^ ( |) P^, + Zg/Zp ^ sin(a)5rf„ (5.26) 

(PdsUd = cos^ {^) Pds - Zs/Zp^smia)Sds, (5.27) 

(5d.)even = COS^ ( |) Srf, - Zp/Zg ^ sin(a)Prf„ (5.28) 

(SdsUd = sm^[^)Sds + Zp/Zs^smia)Pds, (5.29) 



where terms of 0(a) have been neglected. To discuss contributions to the correlation 
function from which — > 7r+ matrix elements are extracted we first assume that the 
pion and kaon source fields have been chosen with definite physical parity, i.e. up to 
cutoff effects they generate excitations with only pion or kaon quantum numbers and 
thus of odd physical parity. It is then clear that the odd parity parts of the renormalised 
operator can only generate 0(a) contributions to the correlation function. Neglecting 
these we are thus left with contributions of the parity even counterterms, i.e. only the 
combination, 

^4(Psd)e.en + ^4(Ssd)e.en^ (5-30) 

contributes at 0(1). Moreover, as can be inferred from eqs. (|5.26I5.28|) . the parity even 
operators are proportional to each other 

iPds)e^en (5'ds)even > (5-31) 

up to terms of 0(a). This means that there is only a linear combination of counterterms 
which contribute at 0(1), which can be written as 

|a4 - Zp/Zg?^^^^Ac|} (P..)_ ^ 4 (P..)_ . (5.32) 

Moreover, it is not necessary to determine the parity even part of Pgd, as the corre- 
sponding parity odd part does only contribute at 0(a). It is therefore enough to insert 
for the renormalised operator 

(Or)±a+av {OvA+AV + dpPsd} ■ (5.33) 

In other words: for insertions into the K ^ n correlation function with meson sources 
of definite parity, the replacement ()5.33|) yields equivalent results up to 0(a) effects. 
The remaining coefficient dp is then determined non-perturbatively by the condition 
(|5.5|) on the hadronic matrix element. 

Finally, by the same arguments as in the preceding subsection we may relax the 
requirement that the meson sources have a definite parity. This is again due to the 
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expected energy gap between scalar and pseudoscalar channels, which implies an ex- 
ponential suppression of the wrong parity contributions. The sources with a definite 
parity therefore just provide an additional suppression by a power of a of these terms, 
which is not necessary for the extraction of the matrix element. 

In conclusion, in the partially twisted case with /3 = 0, the two linearly diver- 
gent counterterms can be subtracted non-perturbatively by imposing two independent 
but otherwise arbitrary conditions on correlation functions involving the four-quark 
operator. Once this has been achieved it turns out that wrong parity contributions 
to the correlation function are of 0(a) provided that the mesonic source fields have 
odd parity. Otherwise they are in any case exponentially suppressed in the asymptotic 
regime where the matrix element is extracted. A remaining finite subtraction constant 
is determined in the usual way by the matrix element itself from eq. (|5.5|) . 

5.6 — > vr'^ matrix elements 

As stated earlier, the breaking of isospin symmetry in twisted mass lattice QCD ob- 
scures the decomposition of four-quark operators into A/ =1/2 and A/ = 3/2 pieces. 
In order to obtain both effective couplings, and §27, one may therefore opt to com- 
pute, besides the —>■ tt"*" transition, also — > vr'^ matrix elements. 

The strategy used for — > vr"*" matrix elements can be taken over virtually un- 
changed. However, a technical complication is generated by the fact that disconnected 
diagrams are now possible, as the self-contraction of the pion source does not vanish by 
flavour symmetry. Moreover, one may be worried about new power divergences, which 
may arise in the pion source field. For instance, the field 



suffers from a quadratic divergence oc ^i/a?. However, we note that this can be easily 
bypassed with the choice of the axial current. 



which, due to the 0(4) vector structure, does not suffer from additive renormalisation. 
Hence we conclude that, beyond possible practical problems relating to the quality of 
the signal in the numerical simulations, the — tt^ matrix elements do not present 
any further theoretical challenge. 

5.7 A comment on K — > tttt transitions 

In principle, there is no obstacle to applying the strategies described here to the physical 
K TTTT amplitudes. Parity is changed in this transition, so that only the parity odd 
operator contributes to the matrix element. This operator has nice renormalisation 



(5.34) 



(5.35) 
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properties already with standard Wilson quarks (|2.28j) . so that the potential gain by 
using twisted mass QCD mainly consists in the elimination of unphysical zero modes. 

A direct mapping of the operator Oy^_|_^y in the twisted basis to the parity odd 
operator in the standard basis is obtained for twist angles a = —(3 = 7r/2 [cf. eq. H3.28() ]. 
Compared to the case of /? = 7r/2 this implies that the sign of and /ic are reversed. 
Alternatively, one may again set /3 = 0, and project on the desired matrix element by 
choosing interpolating fields for the two-pion and kaon states with the correct physical 
parity. However, a practical problem may arise from isospin symmetry breaking in 
lattice regularised tmQCD: for instance, the lattice symmetries do not distinguish the 
/ = two-pion state from the state with a single neutral pion. As the energy of a single 
pion is lower than the energy of the two-pion state, one needs to carry out a multi- 
state analysis which renders this approach more complicated. Nevertheless, the good 
news remains that unphysical zero modes can be eliminated whilst the renormalisation 
properties are no worse than in standard lattice QCD with Wilson type quarks. 

6 Conclusions 

Twisted mass lattice QCD provides a viable framework for a determination of and 
927) by allowing for comparatively cheap numerical computations oi K ^ tt matrix 
elements at light meson masses. In order to facilitate future practical implementations 
of this strategy, we summarise our results in the form of a recipe: 

(i) Choose the number of dynamical quarks to be 0, 2,3 or 4, with mass degenerate 
up and down quarks; the recommended twist angles are then (a,/3) = (7r/2,7r/2) 
or = (7r/2,0). The latter choice avoids a complex quark determinant if 
more than 2 quarks are dynamical. 

(ii) Translate the twist angles to bare mass parameters: (a,/?) = (vr/2,7r/2) means 
that all standard mass parameters are tuned to the critical mass mcr, which can 
be obtained from the PCAC relation as usual. The physical quark masses are 
then determined by the twisted mass parameters fii, fig, fJ'c which are fixed by 
matching appropriate experimental quantities. 

The second option, (7r/2,0), corresponds to = ^''c = 0. For given mass pa- 
rameters ms,mc, one then needs to determine m; so that the axial current in the 
light quark sector is conserved. A simplification occurs with quenched strange 
and charm quarks, where this is equivalent to setting mi = rricr, independently of 
ms,mc. The physical quark masses are determined in terms of ^urug and rric. 

(iii) Choose your preferred mesonic source fields for the kaon and pion; these need 
not be of definite physical parity, however, they should have a sufficiently large 
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overlap with the physical kaon or pion state. In the case of full twist, (7r/2,7r/2), 
there are two ways to proceed: 

• if 0(a) improvement of the massless action and the axial current is imple- 
mented via the counterterms oc Cgw and oc ca respectively, the tmQCD action 
can be considered 0(a) improved and the linear divergence in the four-quark 
operator need not be subtracted. Just compute the bare 3-point function and 
the pion and kaon propagators, and go to the asymptotic regime to obtain 
the matrix element up to 0(a) effects. Requiring the matrix element to van- 
ish at zero meson momenta fixes the coefficient of the finite counterterm. 

• if 0(a) improvement is not implemented, one needs to first determine Cp 
e.g. by imposing the parity restoration condition 1)5. 13(1 . Using the resulting 
subtracted operator one may proceed exactly as above. 

In the case of partial twist, {it/2,0), one first subtracts the linear divergences 
by imposing two arbitrary renormalisation conditions such as (|5.2ni5.21|) . Then 
one proceeds as above, using the subtracted operator. To fix the finite countert- 
erm through the K ^ tt matrix element at vanishing meson momenta, one may 
use the substitution ()5.33|) which avoids determining the parity even part of the 
counterterm. 

(iv) Repeat the calculation for a sequence of lattice spacings, keeping the physical 
scales fixed, and extrapolate to the continuum limit. Note that there exist various 
alternative strategies to compute the same matrix element. By varying the choice 
of twist angles, the way in which the linear divergences are subtracted, or the 
choice of mesonic source fields, one should obtain always the same matrix element 
up to 0(a) effects. Thus, one may hope to improve the control of the continuum 
limit by pursuing a couple of alternative strategies in parallel. 

(v) The continuum results obtained for various meson momenta and energies should 
be well described by leading order chiral perturbation theory, provided the meson 
masses and momenta are small enough. If so, the — > vr"^ matrix element 
allows to determine the linear combination 1^27 + 98) while the matrix element 
for neutral mesons yields §27 — gs (cf. subsect. 2.3). If the AI = 1/2 rule can 
be explained within QCD, one would expect the lattice results to reproduce the 
hierarchy gs » 327 obtained for the phenomenological estimates 1)2. 18() . 

The strategy described above offers major improvements over previous attempts to 
compute K ^ TT transitions with Wilson-type quarks. Nevertheless, it is conceivable 
that further theoretical developments, along the lines of refs. I4 2l35j . may lead to in- 
teresting variations and/or alternatives to the current work. 



32 



Acknowledgments 



We thank H. Wittig for a critical reading of the manuscript, and A. Donini, R. Frezzotti 
and G.C. Rossi for discussions. The hospitality of the Theory Division of CERN (CP.), 
the Physics Department at the University of Rome "Tor Vergata" (CP.), and the 
IFT at the Autonomous University of Madrid (CP., A.V.) is gratefully acknowledged. 
CP. has been partially supported through the European Community's Human Potential 
Programme under contract HPRN-CT-2000-00145, Hadrons/Lattice QCD. 



References 

[1] W.A. Bardeen, A.J. Buras and J.M. Gerard, Phys. Lett. B 192 (1987) 138; 

S. Bertolini, J.O. Eeg, M. Fabbrichesi and E.I. Lashin, Nucl. Phys. B 514 (1998) 
63, and references therein; 

J. Bijnens and J. Prades, JHEP 9901 (1999) 023; 

T. Hambye, S. Peris and E. de Rafael, JHEP 0305 (2003) 027. 

[2] D. Becirevic et al. [SPQcdR Collaboration], Nucl. Phys. Proc. Suppl. 119 (2003) 
359. 

[3] D. Pekurovski and G. Kilcup, Phys. Rev. D 64 (2001) 074502; 
T. Bhattacharya et al., Nucl. Phys. Proc. Suppl. 129 (2004) 257. 

[4] T. Blum et al. [RBC Collaboration], Phys. Rev. D 68 (2003) 114506; 
A. Soni, Pramana 62 (2004) 415; 

J.I. Noaki et al. [CP-PACS Collaboration], Phys. Rev. D 68 (2003) 014501. 

[5] D. Becirevic, Nucl. Phys. Proc. Suppl. 129 (2004) 34. 

[6] L. Maiani and M. Testa Phys. Lett. B 245 (1990) 585. 

[7] L. Lellouch and M. Liischer, Commun. Math. Phys. 219 (2001) 31. 

[8] C. J. Lin, G. Martinelh, C. T. Sachrajda and M. Testa, Nucl. Phys. B 619 (2001) 
467. 

[9] C. Kim and N.H. Christ, Nucl. Phys. Proc. Suppl. 119 (2003) 365. 

[10] C. W. Bernard, T. Draper, A. Soni, H. D. Politzer and M. B. Wise, Phys. Rev. D 
32 (1985) 2343. 

[11] C. W. Bernard, T. Draper, G. Hockney and A. Soni, Nucl. Phys. Proc. Suppl. 4 
(1988) 483. 

[12] C. Dawson et al., Nucl. Phys. B 514 (1998) 313. 



33 



[13] J. Laiho and A. Soni, Phys. Rev. D 65 (2002) 114020; hep-lat/0306035. 

[14] C. J. Lin, G. Martinelli, E. Pallante, C. T. Sachrajda and G. Villadoro, Nucl. Phys. 
B 650 (2003) 301. 

[15] S. Capitani and L. Giusti, Phys. Rev. D 62 (2000) 114506. 

[16] S. Capitani and L. Giusti, Phys. Rev. D 64 (2001) 014506. 

[17] S. Aoki and Y. Kuramashi, Phys. Rev. D 63 (2001) 054504. 

[18] M. Bochicchio, L. Maiani, G. Martinelh, G. C. Rossi and M. Testa, Nucl. Phys. B 
262 (1985) 331. 

[19] L. Maiani, G. Martinelh, G. C. Rossi and M. Testa, Nucl. Phys. B 289 (1987) 505. 
[20] P. Hernandez and M. Laine, JHEP 0301 (2003) 063. 

[21] L. Giusti, P. Hernandez, M. Laine, P. Wcisz and H. Wittig, hcp-lat/0407007. 

[22] C. Pcna, S. Sint and A. Vladikas, Nucl. Phys. Proc. Suppl. 119 (2003) 368; Nucl. 
Phys. Proc. Suppl. 129 (2004) 263. 

[23] A. Pich, B. Guberina and E. de Rafael, Nucl. Phys. B 277 (1986) 197; A. Pich 
and E. de Rafael, Nucl. Phys. B 358 (1991) 311. 

[24] J. Gasser and H. Leutwyler, Annals Phys. 158 (1984) 142; Nucl. Phys. B 250 
(1985) 465. 

[25] R.J. Crewther, Nucl. Phys. B 264 (1986) 277. 

[26] J. Gasser and U.G. Meissner, Phys. Lett. B 258 (1991) 219. 

[27] C. Itzykson and J. B. Zuber, "Quantum Field Theory," McGraw-Hill (1980), New 
York, U.S.A. 

[28] M. B. Gavela, L. Maiani, S. Petrarca, G. Martinelh and O. Pene, Phys. Lett. B 
211 (1988) 139. 

[29] A. Donini, V. Gimenez, G. Martinelh, M. Talevi and A. Vladikas, Eur. Phys. J. C 
10 (1999) 121. 

[30] R. Prezzotti, P. A. Grassi, S. Sint and P. Weisz, Nucl. Phys. Proc. Suppl. 83 (2000) 
941. 

[31] R. Frezzotti, P. A. Grassi, S. Sint and P. Weisz [ALPHA Collaboration], JHEP 
0108 (2001) 058. 

34 



[32] R. Prezzotti and S. Sint, Nucl. Phys. Proc. Suppl. 106 (2002) 814. 

[33] M. Liischer, S. Sint, R. Sommer and P. Weisz, Nucl. Phys. B 478 (1996) 365. 

[34] T. Kaneko et al. [CP-PACS Collaboration], Nucl. Phys. Proc. Suppl. 129 (2004) 
188. 

[35] R. Prezzotti and G. C. Rossi, Nucl. Phys. Proc. Suppl. 128 (2004) 193. 

[36] T. Bhattacharya, R. Gupta, W. J. Lee and S. R. Sharpe, Nucl. Phys. Proc. Suppl. 
83 (2000) 902. 

[37] T. Bhattacharya, R. Gupta, W. J. Lee, S. R. Sharpe and J. M. S. Wu, Nucl. Phys. 
Proc. Suppl. 129 (2004) 441. 

[38] R. Prezzotti, S. Sint and P. Weisz [ALPHA Collaboration], JHEP 0107 (2001) 
048. 

[39] G. Martinelh, C. Pittori, C. T. Sachrajda, M. Testa and A. Vladikas, Nucl. Phys. 
B 445 (1995) 81. 

[40] K. Jansen et al, Phys. Lett. B 372 (1996) 275. 

[41] M. Guagnelli, J. Heitger, C. Pena, S. Sint and A. Vladikas [ALPHA Collaboration], 
Nucl. Phys. Proc. Suppl. 106 (2002) 320; Nucl. Phys. Proc. Suppl. 119 (2003) 436. 

[42] R. Prezzotti and G. C. Rossi, JHEP 0408 (2004) 007; hep-lat/0407002. 



35 



